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ABSTRACT 



Recently, Owen(1993) and Owen and Ledlow (1994) have shown that the 
dividing line between Fan aroff- Riley class I and class II Radio Sources is very sharp 
when the sources are represented as points in the radio-optical luminosity plane. 
It is shown that if one accepts the propositions that the sources in the vicinity 
of this dividing line are characterized by (1) transonic Mach numbers and (2) 
mildly relativistic velocities (consistent with deceleration of an initially moderately 
relativistic or ultrarelativistic jet), then the slope of the dividing line is readily 
obtained using simple physical arguments and established empirical relationships 
between the X-Ray luminosities, core radii, velocity dispersions and absolute 
magnitudes of elliptical galaxies. The intercept of the dividing line depends upon 
parameters which are known perhaps to within factors of order unity and agrees 
with the data to within an order of magnitude. ROSAT observations of elliptical 
galaxies will be important in constraining the central pressures and X-Ray core 
radii of radio ellipticals. Knowledge of these two parameters will assist in a more 
detailed assessment of the physics which is proposed here as being relevant. High 
resolution observations of jets within a kpc of the radio core will also be useful for 
determining the spreading rates of jets in this region. 
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1. Introduction 

Fanaroff and Riley (1974) showed that Radio Galaxies exhibit a change in morphology from 
edge-darkened to edge-brightened at a monochromatic power corresponding to approximately 
10 24 - 5 W Hz -1 at a rest frame frequency of 1.4 GHz. This change occurs over about a decade of 
radio luminosity. It now appears to be well accepted that the mode of energy transport in the 
two types of radio source are different with class-I sources containing turbulent transonic jets 
and class II sources containing jets which are super/hypersonic. 

Owen (1993) and Owen and Ledlow (1994) have recently published new observational data 
which relates to the Fanaroff- Riley classification in a fundamental way. Their survey of cluster 
radio sources shows that, when radio sources are plotted as points in the Radio luminosity - 
Optical luminosity plane, the Fanaroff- Riley break becomes very sharp with the break radio 
power approximately proportional to the optical luminosity squared. (Their data are plotted 
in figure 5). The aim of this paper is to show that this dividing line can be understood on 
the basis of some fairly simple physics together with the assumption that class I jets are at 
least moderately relativistic on the pc scale combined with established empirical relationships 
between core radius, velocity dispersion and absolute magnitude. The assumption of relativistic 
pc-scale flow in FRI sources is supported by a number of recent papers (e.g. Laing, 1993; 
1994; Bicknell, 1994a, b: Parma et al 1994) which show, on several grounds, that it is quite 
feasible that the jets in class I sources are initially relativistic (or at least mildly relativistic) 
decelerating to subrelativistic flow in the first few kiloparsecs. 

De Young (1993) has previously given some consideration to the physical ramifications of 
the Owen-Ledlow data and has used a combination of empirical and theoretical results on the 
transition of jets to turbulence to show that deceleration of initially supersonic jets within the 
core of the parent galaxy is possible. This result is interesting since the deduced transition 
length is comparable to the "gap" lengths observed in some class I sources near the FRI/II 
break. De Young, however, did not consider the implications of deceleration for the distribution 
of radio galaxies in the radio-optical plane and it is this distribution which is the subject of 
this paper. A significant difference between the De Young theory and that developed here, is 
that De Young deduces a much lower velocity (~ the ISM sound speed) for class I jets when 
they become transonic. This is related to the assumed rapid transition of the jet density to 
that of the background ISM. This in turn requires a large ISM number density (n ~ 10 cm -3 ). 
In the theory developed in the present paper the velocities of borderline FRI/II jets are mildly 
relativistic close to the core. 

In the following section the energy flux in a relativistic jet and the ramifications of 
decelerating an initially relativistic jet to a transonic Mach number are considered. 

2. Jet Power 

In this section the basis for the expression for the FRI/II dividing line, the expression 
for the jet energy flux, is discussed. It is also shown that when an initially relativistic jet 
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decelerates, a more or less unique velocity (~ 0.6 — 0.7 c) is obtained. This has the effect, when 
the results of later sections are taken into account, of making the jet energy flux almost entirely 
dependent upon environmental parameters. 

The energy flux of a jet whose internal energy is dominated by relativistic particles is 



where pj et , rj et , 7j et and fj et refer to the jet pressure, radius, Lorentz factor and /3j et = fj e t/ c - 
The parameter 1Z = pc 2 /(e + p) is the ratio of cold matter energy density to enthalpy. [See 
Bicknell (1994b) for the origin of this form of the energy flux and the significance of 1Z]. Thus, 



where the logarithm is base 10. In this equation, as in all similar subsequent equations, terms 
are grouped in three categories, those which involve fixed parameters or constants, those which 
contain parameters which are expected to vary minimally and which may be constrained to 
say, within a factor of order unity, and parameters which are known to vary significantly with 
absolute magnitude. 

The transition point between transonic and supersonic flow is taken to be Ai ~ 2. As 
stated above and in § 1 it is also assumed that the jet flow on the pc scale, (at least for jets near 
the transition between class I and class II) is moderately relativistic (Lorentz factors ;> 2 and 
a not overwhelming contribution to the jet energy density from cold matter) and that the jet 
decelerates between the parsec and kiloparsec scales. Deceleration of a jet in this region may 
occur for two reasons, deceleration due to turbulent interaction with a slow moving, confining, 
nuclear wind {e.g. Smith, 1993) or injection of mass through stellar mass-loss along the path of 
the jet. Both Phinney (1983) and Komissarov (1994) have shown that the latter effect can be 
important when the jet energy flux <J 10 43 erg s _1 . 

In Bicknell (1994b) the velocity of an initially relativistic jet which decelerates to a Mach 
number ~ 2 was shown to be ~ 0.7 c. In that paper it was assumed that the jet inertia is 
dominated by the inertia of the ultrarelativistic particles, that is, there is little cold matter. It 
was also assumed that the jet is initially free. However the calculation is easily extended to 
confined jets in which cold matter may or may not be present. To allow for cold matter (e.g. 
non-relativistic protons) the relativistic enthalpy, w can be expressed in the form: 



where p is the density of cold matter and the internal energy density, e, and pressure, p, are 
assumed to be dominated by relativistic particles. Neglecting the effect of buoyancy on the jet, 
so that the momentum flux is conserved, the equations for momentum and energy conservation 
become, respectively [see Bicknell (1994b) for the basis of these]: 





(2-1) 



w = pc 2 + e + p = pc 2 + 4p 




[(72 - 1)^2 + 72] 72 fh = [(71 - + 7i] 7i fh 
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where the subscripts 1 and 2 refer to the parsec scale and kiloparsec scale respectively, A 
is the jet cross-sectional area and Ap is the excess of jet pressure over the ambient value. 
It is assumed that the kpc scale jet is in pressure equilibrium with the ISM. The quantity 
(pi Ai)/(p2 A2) may be eliminated from these equations to give the following relationship 
between 1Z = pc 2 /4p and (3: 

U2 ~ 1-72(1-/32^) (2 " 2) 

where 

7i(7i-l)/3i^i + 7i 2 /3i 
1 (l + ^ 1 ) 7l 2 /3 1 2 + A Pl /4p 1 

and the relativistic Mach number (AA = (2 + 37?.) 1 / 2 7/3) follows from this. 

The relationship between Mach number and velocity is shown in figures 1 and 2 for both 
confined and unconfined jets, for a range of initial ratios of cold matter energy density to 
enthalpy, and for different initial Lorentz factors ranging from 1.5 to 10. In figures 3 and 4, 
the corresponding relationship between 7^2 an( i velocity is also plotted. Clearly, although the 
details of the curves change, the result is approximately the same: When an initially relativistic 
jet decelerates to a Mach number of 2 and cold matter does not greatly dominate the jet 
inertia, its velocity lies approximately between 0.6 and 0.7 c. Parma et al. (1994) have shown, 
for a sample of Bologna sources, that such velocities are consistent with surface brightness 
asymmetries at the bases of class I jets. 

The Mach number-velocity curves are accompanied by a somewhat curious behavior in the 
parameter 1Z. When this parameter is initially high, deceleration of the jet (that is, decreasing 
P2) implies that 1Z decreases. This is probably related to the fact that the sound speed in the 
high 1Z jets is lower, the Mach number is therefore higher for the same velocity and if such a 
jet decelerates, it dissipates more kinetic energy into internal energy, lowering the value of 1Z. 
This is not to say that a high Mach number jet will decelerate in such a fashion. This is merely 
the result if it decelerates. 

Note also that the curves for the initially unconfined, high 1Z jets, start (at the right hand 
extremum of the curve where /?2 = fi\) with a higher value of 7^2 than the pc-scale value. This 
is probably related to the fact that such a jet would initially accelerate due to the influence 
of the high pressure and that if it is to decelerate to its initial velocity at a point where it is 
confined, then entrainment must occur. These indirect explanations are due to the fact that 
the description of entrainment here is incomplete and one is relying on inferences based upon 
energy and momentum conservation. 



3. Jet Pressure and Radius at the Transition Point 



3.1. Jet Radius 



In some well-observed galaxies [e.g. IC4296, Killeen et al (1986) and NGC 315, Venturi et 
al., (1993)] in which the jets exhibit class I morphology, the transition to turbulent flow appears 
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to occur not far outside the optical core of the parent galaxy. There are good physical reasons 
for this: This is where the density of the ISM is highest and where the pressure gradient is 
flattest (neglecting for the present, the possible transition region between a nuclear wind and 
the static ISM) and the jet is not driven by the pressure gradient of the ISM. Thus, in using 
equation (2-1) the jet pressure p t at the transition point between supersonic and transonic flow 
is normalized by the core pressure. 

There are several factors which may affect the jet radius in the region before the transition 
to fully turbulent, transonic flow: 

1. The jet expands as described by the equations of mass-conserving supersonic, laminar 
flow. 

2. The jet, although supersonic may expand due to turbulence, albeit at a lower rate than a 
transonic or subsonic jet. This may occur if there is an atmosphere in the core region of 
the galaxy, confining the jet. Such an atmosphere may be either static or in the form of a 
slowly moving wind. 

3. The initial 100 pes or so of the jet may represent a transition region to fully turbulent 
flow. This is the idea developed by De Young (1993) and, as stated in the introduction, 
it is of interest that the transition lengths that De Young calculates are typically of order 
the "gap"-lengths of some borderline class I/class II jets. 

4. Mass injection from the stars in the core of the galaxy as suggested by Phinney (1983) 
and Komissarov (1994). 

A small number of jets have been observed with enough resolution in the region ~ 100 pc 
from the nucleus to show that some jets expand slowly at first and then suddenly start 
to expand more rapidly. For example, both of the IC4296 jets initially expand at a rate 
d$/dQ « 0.2, (Killeen, Bicknell k Ekers, 1986) the NGC315 main jet initially expands with 
d^/dQ ~ 0.1 (Bridle, 1982) . The observations of such jets suggest that something other than 
mass conserving supersonic laminar flow is involved in the initial region. However, as outlined 
above, the physics of this expansion clearly involves a number of different possibilities which 
are not explored here. In the following, a fiducial value of rj et /r t = 0.1 is adopted where r\ is 
the distance from the core at which the transition to turbulent flow occurs and rj et is the jet 
radius at that point. A value for this parameter ~ 0.1 is also justified by other observations: 
Bridle's (1985) compilation of jet mean spreading rates (< d^/dQ >) as a function of core 
power shows a dramatic change at the class I-class II boundary from approximately 0.3 to less 
than 0.1. If we regard the transition spreading rate (~ 0.1) as typical of supersonic jets with 
Mach numbers just above transonic, then this would imply rj et /r t ~ 0.05. 

3.2. Jet Pressure 

The core pressure of the ISM may be related to the X-Ray luminosity of the galaxy as 
follows. Let us represent the number density, ra(r), of the hot galactic atmosphere by the 
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following empirical relation, frequently used in analyses of X-ray data: 



n(r) = n c 



r 2 



"/3a 



(3-1) 



where n c is the central number density, r c is the core radius of the X-ray emitting gas. 

Contrary to some earlier uses of this relation, this does not represent an isothermal 
distribution of gas in an isothermal sphere (see Killeen and Bicknell, 1988 for a discussion 
of this point). However, the spherically symmetric equations for an isothermal hydrostatic 
atmosphere and the gravitating matter distribution, namely, 

1 d Pg = ( 3 _2) 
p g dr dr 

1 d ( 2 \ d<f) 



^Tr^) = ~Tr (3 " 3) 



d 2 c/) 2 d(j> 
dr 2 r dr 



4irG Pm (3-4) 



may be used to derive the details of the matter density and velocity dispersion implied by 
equation (3-1). (In the above the subscript 'g' refers to the gas and 'm' refers to the gravitating 
matter; (f> is the gravitational potential, ay is the radial component of velocity dispersion and 
G the constant of gravitation.) 

The matter density follows from using equation (3-2) to derive the potential gradient 
followed by insertion into Poisson's equation (3-4). The stellar hydrodynamic equation, 
(3-3) then yields the velocity dispersion. The gas is assumed to be isothermal with 
p g = (pkT)l(jim v ). The results are: 



°lo 3 + r 2 /r 2 c 



2irGr 2 (l + r 2 /r 2 ) 2 
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2 + r 2 /r 2 c 



00 3 + r 2 /r 2 



at 



Gr c (1 + ^2/^2)3/2 

kT 



The parameter is the radial velocity dispersion at infinity and if the velocity dispersion is 
isotropic, this is the observed value at large radii; S(i?) is the surface density as a function of 
projected radius R; the last relation between /3 at and the velocity dispersion is a consequence 
of ensuring that the stellar kinetic energy density vanishes at large radii. 

While the above "law" for the gas distribution is empirical, the implied density profile 
does have some appealing features consistent with a distribution dominated by luminous 
matter within, say an effective radius and then merging into an isothermal sphere, dark matter 
dominated distribution at large radii. First note that the surface density drops to 0.5303 
of its central value at R = r c so that the parameter r c can be equated to the optical core 
radius. Second, the central velocity dispersion is smaller in the center compared to the outside, 
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mimicking to some extent, the behavior inferred for some ellipticals by Fall (1987). In the 
above distribution = 1.5 Cg where <jq is the central velocity dispersion. Fall's compilation of 
data on early type galaxies shows that (with a wide scatter) the mean circular velocity, v c , is 
approximately 2.1 &q. With an isothermal distribution at large radii, v c = \/2g 00 and therefore, 
(7^, = 2.2 Gq. Hence the data reviewed b Fall and the above model are not totally in agreement. 
However, the discrepancy is not too bad and this model is a reasonable one with which to 
estimate the central pressures of X-Ray emitting atmospheres. High resolution ROSAT data 
can be expected to yield further interesting information on the distribution of the ISM within 
ellipticals. 

The temperature of the X-ray emitting gas is given by: 



kT 



fj, m, p g c 



2.2 



Ait A.t 

the latter following from Fall's (1987) relation between circular velocity and central velocity 
dispersion. 

The X-ray luminosity in the Einstein 0.2-4 keV band implied by equation (3-1) is 

Lx = 4vr (^j n e 2 iC r c 3 A E (T)/ x (Ait,ri/r c ) 
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where 
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n t is the total ion density, n e is the electron density Ae is the cooling function integrated over 
the Einstein band, T is the temperature of the (assumed isothermal) X-Ray emitting gas and 
r\ is a cutoff radius. The values of n t /n e and Ae presented in table 1, based upon calculations 
using MAPPINGSII (Sutherland and Dopita, 1993), were kindly provided by R. Sutherland. 
Obviously there is not a great change in Ae over the range of temperatures of interest and the 
value Ae ~ 1.2 X 10~ 23 appropriate to T = 10 7 K is used. 

Collecting the above results, the central pressure of the X-Ray emitting atmosphere 
[p c = (n t + n e ) kT] can be expressed as: 

n t /n e 



fogPc 
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log(2.2/xmp) - 2 lo S77 
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(3-7) 



Typically, models of X-ray emitting atmospheres imply a mean value for the parameter 
Ait ~ 0-75 (e.g. Thomas et al., 1986). For /3 at < 0.75 the integral Ix diverges (logarithmically 
for Ait = 0.75). Values of ix are given in table 2 for various values of /3 at and r\jr c . 



Substituting equation (3-7) for the pressure in equation (2.) for the energy flux, gives: 
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4. Summary of Empirical Relationships between X-Ray Luminosity, Velocity 

Dispersion and Absolute Magnitude 



Clearly, the expression for the energy flux derived in the previous section depends upon 
the X-ray luminosity, the velocity dispersion and the X-ray core radius, all of which are 
empirically related to the optical luminosity. Here the empirical relationships between the 
relevant parameters are summarized. All of the following relationships have been transformed 
to a uniform value of the Hubble constant, Hq = 75 kms -1 Mpc -1 . 



4.1. Relationship between X-Ray and Optical Luminosities 



Donnelly et al (1990) have reconsidered the correlation between X-Ray and optical 
luminosities for ellipticals using revised distances. Their correlation can be expressed in the 
form: 

log i x = 22.3 - 0.872 M B (4-1) 



4.2. Relationship between Core-radius and Absolute Magnitude 



This has been summarized by Kormendy (1987). The best fit straight line to the data on 
optical core radii is given by: 

log r c = 11.7 - 0.436 M B (4-2) 
As indicated in the previous section, I take the optical and X-ray core radii to be identical. 



4.3. Relationship between Velocity Dispersion and Absolute Magnitude - the 

Faber-Jackson Relation 



A fairly recent determination of the Faber-Jackson relation (Terlevich et al., 1981) yields, 
for H = 75, 

log a = 5.412 - 0.0959 M B (4-3) 
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where the velocity dispersion, <7, is in units of cms 



5. The Class I/Class II Dividing Line 



5.1. The Central Pressure 



Insertion of the above correlations into the expression (3-7) for the central pressure gives: 

logp c = -8.2 + 0.03 M B - i log 

showing a minor dependence on absolute magnitude. [Such an effect had been remarked 
upon by Thomas et al. (1986).] For a typical Mb = —21 this gives a central pressure 
p c f» 2 X 10~ 9 dy cm -2 . 




5.2. The Jet Energy Flux 



Combining equations (3-8) and the equations relating the X-ray and optical parameters in 
the previous section, one obtains for the jet energy flux: 
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Inserting numerical values for physical constants and the approximately determined parameters 
discussed earlier, 
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(5-1) 



Inserting a typical blue magnitude of -21 for a radio elliptical, taking /3j et = 0.6 and adopting 
fiducial values for the other parameters, one obtains ~ 2 X 10 42 ergs s _1 . To order of 
magnitude this is the energy flux that is often assumed for borderline class I/class II radio jets. 



5.3. The Radio Power 



The final connection required is that between the jet energy flux and the radio power 
at 1.4 GHz. Taking k to be ratio of the non-thermal luminosity of one side of the source to 
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the corresponding jet energy flux we have for the total synchrotron plus inverse Compton 
luminosity (Ls and Lie respectively): 

L s + Lie = 2 k F E 

Moreover, the monochromatic power at frequency vq is given in terms of the synchrotron 
luminosity and the upper cutoff frequency z/ u and spectral index, a, by: 



L V0 = v^Cs (%,<*) L S 
where C s (^,a) = 

is tabulated in table 3. In view of this table, a fiducial value of C = 0.1 is adopted. 
Hence, 



log 
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and for v = 1.4 X 10 9 GHz, 
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Combining this with the expression for the energy flux (equation (5-1) and utilizing the 
mean color difference for ellipticals [< B — R >= 1.9 from tabulated data of Persson et al, 
1990)] gives: 
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(5-3) 



Most of the terms in the expression for log A are arguably known to within factors of order 
unity - save for the ratio k of lobe luminosity to jet energy flux. Following Bicknell (1986) this 
may be constrained in the following way: The energy budget for a lobe is given by 



where El is the lobe energy, P the pressure, V the volume, 
is the lobe luminosity and the radiative time scale 
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where / is the fraction of the lobe energy due to electrons and/or positrons and c\2 is a 
Pacholczyck (1970) synchrotron parameter. This parameter is mainly dependent on the 
upper cutoff frequency, decreasing as the latter increases. If / ~ 0.5 and B ~ 3 [iG then 
W ~ 2 x 10 8 yr; if B m 10 fiG, W ~ 1 X 10 7 yrs. 

For illustration, consider a particle-dominated plasma, then PdV/dt = 0.25 F E (see 
Bicknell, 1986) and 

dE L _ 3 x 

^r-4 FE -^ EL 

Ift» W then clearly, C = t~^E L -> 0.75F E . If t « t Uad then C « 0.75(t/t Tad )F E . The 
details of this calculation will change if the plasma is say, in equipartition, but not the general 
principles. 

In principle, for a source in which / ~ t ra a, the value of k can be of order unity. Whether k 
can approach the limiting value of 0.75 depends upon whether the radiative time scale can be 
kept reasonably short by acceleration of high energy electrons, preventing the value of c\2 (and 
hence i ra d) from increasing. This could occur if the lobes are turbulent. 

The sources NGC 315 and NGC 6251 were examined in detail in Bicknell (1994b) and 
it is of interest what that analysis implies for the respective values of k. The relevant source 
parameters are summarized in table 5.4. Since the magnetic field in the lobes of both of these 
sources is low (-B m i n ~ 1.5/iG) the inverse Compton luminosity is significant and the total 
non-thermal luminosity is taken to be a factor of 5 greater than the synchrotron luminosity. In 
the table, K est is the value of k estimated from k = 0.75 t/t ra a, with Fe, est the corresponding jet 
energy flux in each jet. (Such estimates are only possible for NGC 6251 since there is no age 
estimate for NGC 315.) F E ("f = 2) and k(j = 2) are respectively the energy fluxes required if 
the core jets have Lorentz factors ~ 2 and P/P m i n ~ 10 on the pc scale (Bicknell, 1994). 

As the table shows the relevant values of k in these two sources are of the order of 
0.05 — 0.1 if their core jets are moderately relativistic. In the case of NGC 6251 a value ;> 0.05 
is supported by the age of the source estimated from spectral steepening. It is important to 
note that both of these sources are large, presumably due to the fact that they reside in poor 
groups. Sources confined in richer environments (as in the Owen-Ledlow cluster survey) would 
be expected to have larger magnetic fields, and a correspondingly smaller radiative time scale, 
t ra d- It should also be noted that Parma et ol. (1994) found that doppler induced brightness 
asymmetries at the bases of class I jets and estimates of jet velocities based upon the energy 
budget are consistent when k ~ 0.05. In the following a fiducial value of 0.1 is adopted for k. 

5.4. Comparison with Owen-Ledlow Data 

Figure 5 shows the Owen-Ledlow data with two theoretical lines defined by equations (5-2) 
and (5-3) overlaid. These lines correspond to values of k equal to 0.1 and 0.75. The other 
parameters used to define this line (via the equation for log A) are summarized in table 5.4.. 
The value of Lie I Ls was taken to be zero on the grounds that a more tightly confined cluster 
source (as in the Owen-Ledlow sample) would have a larger magnetic field. As one can see, 



the theoretical lines reproduce the slope of the division between class I and class II sources 
very well. The line for k = 0.75 which represents the upper limit for this parameter bisects 
the class I and class II data quite well. The line for k = 0.1 (which perhaps represents a 
lower limit on this parameter) is almost an order of magnitude in log Lr below the actual 
dividing line indicated by the data. However, this is acceptable given the uncertainty in some 
of the parameters that contribute to log A which defines the intercept of the line. Perhaps, the 
remarkable point about the theoretical relationship is that for a reasonable range of values of 
k the intercept agrees with the data to within an order of magnitude without forcing of the 
other parameters. Nevertheless, let us examine some of the terms which contribute to log A 
and determine in what direction they may vary. 

The term (pt r t) / (Pc r c) has been set to unity. This is valid for an atmosphere in which 
p ex r~ 2 . However, for the sort of atmosphere we have been considering here, well outside a 
core radius 



If the transition to transonic flow does not take place until about ten core radii then this 
parameter is approximately 1.8 for /3 at = 0.75. 

The value of Ix decreases (increasing log A) if the radial cutoff in the atmosphere is smaller 
or if the effective /3 at is greater and perhaps here one should use an extreme value rather than 
an average value since this would maximize the radio power, holding other parameters constant. 
An increase of /3 at to 1.0 would increase logi^ by approximately 0.2. However, perhaps a 
more important effect is that of cooling on pressure in the central region: For an hydrostatic 
atmosphere it is straightforward to show that 



where r\ is a reference radius, the subscript c refers to the center, W is the dimensionless 
potential — ^/u 2 , a is a reference velocity dispersion (say, the velocity dispersion at large radii) 
and (3 = (fim p a 2 )/(kT) as before. The second factor on the right hand side represents the 
contribution of an isothermal atmosphere and the first represents the effect on the pressure due 
to the departure from isothermality. Clearly, the central pressure increases above the value it 
would have in an isothermal atmosphere if the gas cools towards the center. An increase in 
pressure would increase the energy flux and hence the radio luminosity. 

The other major parameter affecting the estimate of the energy flux is the factor 
1 + 7 _1 (7 — l)TZ)j 2 (3 occurring in equation (2-1) for the energy flux. As the value of (3 
increases from 0.6 to 0.8 this parameter also increases by a factor of 1.8. 

Thus, the effect of any likely change in the jet and ISM parameters is to increase the jet 





(5-5) 
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energy flux, producing better agreement between the theoretical and observed dividing lines. 

6. Discussion 

In this paper I have shown that the division between Fan aroff- Riley Class I and Class II 
sources in the radio-optical plane may be understood in terms of the following: 

• Class I jets are transonic or subsonic on the kpc scale. However, the jets which are 
borderline class I or class II decelerate to transonic near an optical core radius. 

• Class I jets are relativistic on the parsec scale but the jets in sources on the class I/class II 
boundary decelerate to mildly relativistic (~ 0.6 — 0.7 c) when they become transonic, 
consistent with energy and momentum conservation and consistent with constraints 
derived from surface brightness asymmetries at the bases of class I jets. 

• The derived relationship between the central ISM pressures of elliptical galaxies and the 
absolute magnitude as well as the dependences of core radius and velocity dispersion on 
absolute magnitude. 

The theoretical dividing line has, within the errors, exactly the required slope and the range in 
the calculated zero point lies within an order of magnitude of that defined by the data. 

There is little scope to change the slope of the dividing line since the parameters which 
contribute to it are well determined. Hence the fact that this slope agrees so well with that 
observed is a positive feature of the theory. On the other hand the parameters which contribute 
to the intercept of the dividing line are not as well determined. However, a good proportion 
of the preceding analysis has been devoted to a discussion of the relevant parameters and as 
a consequence it is clear which parameters are required to be better determined in order to 
obtain a better comparison between theory and observation. 

The central ISM pressure is one of the most important parameters which affects the 
intercept and it has been shown in § 5 how cooling within the central region of the galaxy 
would increase its value. This, in turn, would increase the energy flux and make it more likely 
that closer agreement between theory and observation will be obtained without unrealistic 
assumptions on the ratio of radio power to jet energy flux. ROSAT observations will no doubt 
constrain the central pressures and the relevant X-Ray core radii much better than at present. 
Similarly, increasing the jet velocity from 0.6 to 0.8 c increases the jet energy flux by almost 
a factor of 2. The spreading rate of jets close to the core also affects the energy flux quite 
markedly and therefore a detailed radio study of the spreading rates of kpc scale jets close 
in this region would be of interest. The other parameter which is not well determined is the 
ratio of non-thermal luminosity to jet energy flux. As discussed in the previous section this 
could well be higher in cluster sources than in sources in weak groups such as NGC 315 and 
NGC 6251. 
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While there is some doubt as to the value of the intercept of the dividing line, there 
seems to be little doubt that the intercept fundamentally depends upon the initial velocities of 
the jets on the kpc scale being at least mildly relativistic. The best way to realize this is to 
consider the factor (1 + 7 _1 (7 — l)lZ)j 2 ft in the expression for the energy flux. For (3 ~ 0.6 
and Ai ~ 2 this term is approximately 1.3. For (3 ~ 0.1 the term is smaller by about an order 
of magnitude and the theoretical dividing line would be well below that observed. Moreover, it 
is the existence of a more or less unique velocity which effectively separates environmental and 
jet parameters in the expression for the energy flux and makes possible such a sharp division 
between class I and class II sources. The slope of the dividing line is due exclusively to the 
dependence of central pressure, core radius and velocity dispersion on absolute magnitude. 

One parameter which has been introduced here in an empirical fashion is the pc-scale to 
kpc-scale spreading rate. As discussed in § 3 this could depend upon a number of pieces of 
physics. Investigation of the mechanism for jet spreading in the region near the core and its 
relationship to the atmosphere of the central object remain interesting problems despite the 
insights of Phinney (1983), De Young (1993) and Komissarov (1994). 
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Tables 



T (°K) 


n t /n e 


n 


A E 




5.01 x 10 6 


0.909 


0.620 


1.267 X 10" 


-23 


1.00 x 10 7 


0.909 


0.620 


1.170 X 10" 


-23 


3.16 x 10 7 


0.909 


0.620 


1.094 x 10" 


-23 



Table 1: The cooling function, Ae, integrated over the Einstein 0.2-4.0 keV band together with 
related parameters. 



ri/r, 





10 


100 


1000 


00 


0.5 


8.53 


98.4 


998 


00 


0.6 


4.58 


24.3 


103 


00 


0.7 


2.60 


7.22 


14.6 


00 


0.75 


2.00 


4.30 


6.60 


00 


0.8 


1.57 


2.73 


3.46 


4.72 


0.9 


1.01 


1.32 


1.40 


1.42 


1.0 


0.686 


0.775 


0.784 


0.785 



Table 2: Values of the integral ix used for evaluating the relationship between central pressure 
and X-Ray luminosity. 





Upper 


cutoff frequency (z/ u ) 


a 


10 GH: 


z 100 GHz 


0.7 


0.17 


0.083 


0.8 


0.13 


0.085 


0.9 


0.082 


0.065 



Table 3: Values of the function, C's which relates monochromatic and total luminosity (for 
v = 1.4 GHz.) 
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Parameter 


NGC 315 


NGC 6251 


min 


o X 1U ergs 


o X 1U ergs s 


r „ 


z x iu erg s 


z x iu erg s 


£>s + C>ic 


, , i n42 orn . c l 

~ 11) CI D 


, , i n42 orn . c l 


W = E^J{L S + Cic) 


~ 3 X 10 8 yr 


2 x 10 9 yr 


t 




£ 1.3 X 10 8 yr b 


K est = 0.75 (i/W) 




^0.07 


-p£J,est 




~ 8 x 10 42 ergs" 1 


F E {l = 2) 


~ 5 X 10 42 erg s -1 


~2x lO^ergs" 1 


«(7 = 2) 


~ 0.1 


~ 0.05 



Table 4: Source parameters and estimates of the parameter k. Data for NGC 315 from 
Willis et al. (1981), for NGC 6251 from Waggett, Warner and Baldwin (1977) and Saunders et al. 
(1981). 



"Estimated from regions 9, F, G and H of the maps in Willis et al., 1981 

b The age is underestimated if reacceleration is occurring in the lobes of this source. 



Parameter Value 

rjet/n 0.1 

(Ptrl)l(pcjl) 1.0 

/3,7 0.8,5/3 

I x 5.0 

k 0.1, 0.75 

C s 0.15 

Lic/Ls 0.0 

Table 5: Parameters used for theoretical ClassI/II dividing line. 
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Figure Captions 

Figure 1. The relationship between Mach number and jet (3 for decelerating confined jets with 
initial Lorentz factors 1.5, 2.5, 5 and 10.0 (respectively solid, short dash, medium-dash and 
long-dash curves) and varying initial amounts of cold matter (1Z = 0, 1, 2 and 5). 

Figure 2 The relationship between Mach number and jet (3 for decelerating initially free jets 
with initial Lorentz factors 1.5, 2.5, 5 and 10.0 (respectively solid, short dash, medium-dash 
and long-dash curves) and varying initial amounts of cold matter (1Z = 0, 1, 2 and 5). 

Figure 3 The parameter 1Z = pc 2 /(e + p) as a function of jet (3 for decelerating confined jets 
with initial Lorentz factors 1.5, 2.5, 5 and 10.0 (respectively solid, short dash, medium-dash 
and long-dash curves) and varying initial amounts of cold matter (1Z = 0, 1, 2 and 5) 

Figure 4 The parameter 1Z = pc 2 /(e + p) as a function of jet (3 for decelerating initially free 
jets with initial Lorentz factors 1.5, 2.5, 5 and 10.0 (respectively solid, short dash, medium-dash 
and long-dash curves) and varying initial amounts of cold matter (1Z = 0, 1, 2 and 5). 

Figure 5 Theoretical expressions for the class I/II dividing line overlaid on the radio-optical 
data of Owen and Ledlow (1994). The solid line corresponds to k = 0.1; the dashed line to 
k = 0.75. 
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